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For any D-dimensional quantum lattice system, the fidehty between two ground state many- 
body wave functions is mapped onto the partition function of a D-dimensional classical statistical 
vertex lattice model with the same lattice geometry. The fidelity per lattice site, analogous to the 
free energy per site, is well defined in the thermodynamic limit and can be used to characterize the 
phase diagram of the model. We explain how to compute the fidelity per site in the context of tensor 
network algorithms, and demonstrate the approach by analyzing the two-dimensional quantum Ising 
model with transverse and parallel magnetic fields. 
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The discoveries of high- Tc superconductors and frac- 
tional quantum Hall liquids have stimulated a surge 
of activities in the study of quantum phase transitions 
(QPTs) [if. The conventional description of QPTs in 
condensed matter physics is in terms of orders and fluc- 
tuations. The Landau- Ginzburg- Wilson paradigm de- 
scribes symmetry-breaking orders quantified by a local 
order parameter, whose non-zero value characterizes a 
symmetry-broken phase. Continuous QPTs beyond the 
Landau-Ginzburg- Wilson paradigm also exist. They are 
described in terms of the so-called topological/quantum 
orders 01 and are relevant to emergent phenomena in 
strongly correlated electron systems, with nonlocal order 
parameters as a salient feature. 

By using concepts of quantum information science, re- 
cently two new approaches to study QPTs have been 
proposed. They focus on properties of ground state wave 
functions of the quantum many-body system, nar nely en- 
tanglement d i, H, H, 0, H, M and fidelity [13, III d , 
and turn out to be very successful at detecting quan- 
tum critical behaviors. In particular, the entanglement 
entropy is exploited to reveal qualitatively different be- 
haviors at and off quantum criticality [1, 0], whereas 
the fidelity, a measure of distinguishability of states in 
the system's Hilbert space, is shown to be able to cap- 
ture drastic changes in quantum ground states when the 
system undergoes a QPT, regardless of what type of in- 
ternal order is present [ll| . Both approaches have been 
shown to be insightful in the context of already well- 
understood systems, but in practice, when applied to a 
generic system, they still rely on our ability to compute 
certain properties of ground state wave functions, which 
is in general a very difficult task. 

On the other hand, significant progress has also been 
made recently in the classical simulation of quantum 
many-body systems by using a tensor network (TN) to 
represent the wave function. Exarnples of TNs include a 
matrix product state (MPS) [13, O, [l3| for systems in 
one spatial dimension and the projected entangled-pair 
state (PEPS) [16] in two and higher spatial dimensions. 



For systems invariant under translations, particularly ef- 
ficient algorithms have been proposed to compute the 
ground state for infinite systems, both in one 17[ and 
two spatial dimensions, as well as for finite systems 
with periodic boundary conditions (PBC) [3. 

The purpose of this Letter is two-fold. We consider 
a system, either infinite or finite with PBC, defined on 
a Z3-dimensional lattice and such that its ground state 
is invariant under translations (20| . First, we show that 
the fidelity between two ground states can be mapped 
onto the partition function of a D-dimensional classi- 
cal statistical vertex lattice model with the same lattice 
geometry. This is achieved by exploiting the fact that 
the two ground states can be represented in terms of a 
TN where all the tensors are copies of one single tensor. 
The fidelity per lattice site, introduced in Ref. [ll|, is 
naturally interpreted as the free energy per site of this 
Z?-dimensional classical statistical vertex lattice model, 
and as such it is well defined in the thermodynamic limit 
(even though the fidelity itself becomes zero). Second, 
we consider the practical computation of the fidelity per 
lattice site, both for finite and infinite systems, within 
the framework of TN al gor ithms for translationally in- 
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As a result, we obtain a 



viable scheme to determine the ground state phase dia- 
gram of a system without prior knowledge of order pa- 
rameters. We demonstrate the approach by analysing the 
two-dimensional quantum Ising model with both trans- 
verse and parallel magnetic fields. First and second or- 
der phase transitions, as well as stable fixed points, are 
clearly identified. 

Generalities. Consider a finite quantum lattice system 
S" in D dimensions described by a Hamiltonian H{X), 
where A is a control parameter 2l|. For two ground 
states |^/'(Ai)) and \ip{X2)) corresponding to two different 
values Ai and A2 of the control parameter A, the ground 
state fidelity F(Ai,A2) = \{ip{X2)\i^{Xi))\ asymptotically 
scales as i^(Ai,A2) ^ d{Xi,X2)^ , with N the number of 
sites in the lattice. Herejd(Ai, A2) is the scaling parame- 
ter, introduced in Ref. [ll| for one-dimensional quantum 
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systems, which characterizes how fast the fidehty goes 
to zero when the thermodynamic Umit is approached. 
Physically, the scaling parameter d{Xi,X2) is the aver- 
aged fidelity per lattice site, 



ln(i(Ai, A2) 



lim 



lnf(Ai,A2) 

N 



(1) 



which is seen to be well defined in the thermodynamic 
limit even if i^(Ai,A2) becomes trivially zero. It satis- 
fies the properties inherited from the fidelity i^(Ai,A2): 
(i) normalization d(A, A) = 1; (ii) symmetry (i(Ai, A2) = 
d(A2, Ai); and (ui) range < (i(Ai, A2) < 1. Additionally, 
in a finite system we can define a finite-size analogue of 
'^('^ii A2) through 



In dAr(Ai, A2) 



lnj^(Ai,A2) 
N 



(2) 



As argued in Ref. the fidelity per lattice site 

rf(Ai, A2) succeeds in capturing nontrivial information in- 
cluding stable and unstable fixed points along renormal- 
ization group fiows. Specifically, suppose the system S 
undergoes a QPT at a transition point Ac. Then d(Ai, A2) 
exhibits singular behaviors when Ai crosses Ac for a fixed 
A2, or A2 crosses Ac for a fixed Ai. That is, a transition 
point Ac is characterized as a pinch point (Ac, Ac) for 
continuous QPTs: the intersection of two singular lines 
Ai = Ac and A2 — Ac on the two-dimensional surface de- 
fined by d(Ai,A2) as a function of Ai and A2. For first 
order QPTs, d{Xi,X2)) becomes discontinuous (as cither 
Ai or A2 crosses a transition point) [22]. 

Mapping onto a D-dimensional classical statistical ver- 
tex lattice model. As it is well known, there is a remark- 
able mapping from a _D-dimensional quantum system to 
an equivalent {D + l)-dimensional classical system with 



imaginary time as an extra dimension [231 |24| . Here, 
we discuss another mapping, one from the ground state 
fidelity i^(Ai,A2) for a _D-dimensional quantum lattice 
model onto the partition function of a D-dimensional 
classical statistical vertex lattice model. This mapping 
implies that we can take advantage of the whole ma- 
chinery of the transfer matrix formulation in statistical 
mechanics. As we discuss below, it also means that we 
can compute the fidelity per lattice site d(X\ , A2) by ex- 
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ploiting the TN algorithms of Refs. 

To establish this mapping, we recall that any state of 
a quantum lattice system may be represented in terms 
of a TN, such as an MPS for one-dimensional systems or 
a PEPS for systems in D > 2 dimensions [3, As 
a concrete example, let us consider a square lattice on a 
torus with N = x Ly sites, where each site, labeled 
by a vector r = {x,y), is represented by a g-dimensional 
Hilbert space V^^ = £«. A PEPS for a state \iP{X)) con- 
sists of a set of tensors A^^, one tensor per lattice site. 
Each tensor is made of complex numbers ^[^^^^ labeled 
by one physical index s and four bond indices a, /?, 7 and 



6 (in a generic case, there will be one bond index for 
each outgoing link of site r) . The physical index s runs 
over a basis of V^^, so that s — 1, • • ■ ,q, whereas each 
bond index takes Q values, with Q some inner dimension 
of bonds in the valence bond picture, which connects 
the tensors in the nearest neighbor sites. In terms of 
the PEPS representation, the ground state fidelity turns 
out to be equivalent to the partition function of a two- 
dimensional classical statistical vertex lattice model, see 
Fig. ll]), with the statistical "weights" 



E 



A 



[r\s 



(Ai), (3) 



where the tilded indices are combined pairs of indices: 
a = {a, a') and so on. By inspecting definitions ^ and 
([2]), one concludes that the logarithm of dN{Xi,X2) is 
formally equivalent to the free energy per site in the two- 
dimensional classical statistical vertex lattice model 
(up to an irrelevant prefactor linear in temperature). 
This argument is valid for any lattice geometry in any 
dimensions [1^. Therefore, the fact that QPTs may be 
detected as singularities in d{Xi, A2) matches the conven- 
tional wisdom that phase transition points are refiected 
as singularities, in the thermodynamic limit, of the free 
energy for classical systems. 

Some remarks are in order. First, the mapping is exact 
both for finite lattices (possibly for a large Q) and infinite 
lattices (infinite Q). Second, for periodic systems that 
are invariant under translations, one can always build a 
TN where all the tensors are the same (often at the cost 
of increasing Q) by using results in [l^i US] and gener- 
alizations thereof. Finally, in practical computations as 
described below, the exact ground state is approximated, 
in a controlled way, by a TN with reasonably small Q. 





FIG. 1: (color online) Diagrammatical representation of sev- 
eral tensor networks. Left: two-dimensional tensor network 
for the ground state fidelity _F(Ai,A2) in a system defined 
on a torus. Right: matrix product operator (MPO) for the 
corresponding one-dimensional transfer matrix r(Ai, A2), and 
matrix product state (MPS) for the left and right eigenvectors 
of T, and \^r), with the largest eigenvalue fi. 

Fidelity per lattice site from tensor network represen- 
tations. From now on we specialize to a ZJ-dimensional 
lattice system that is invariant under translations by one 
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lattice site [20[. We explain how to obtain the fidelity 
per site, both in infinite and finite (but lar ge) systems. 
As a first step, we use the TN algorithms [13, El, El 
to compute a TN representation for the ground states 
|'0(Ai)) and |V'(A2)) in terms of site-independent tensors 
AR(Ai) and A^^{X2), that we use to build the (also site- 
independent) statistical weig hts E^^{Xi,X2). We notice 
that all these tensors depend on the lattice size A^. 

The fidelity -F(Ai, A2), regarded as the partition func- 
tion of a _D-dimensional classical statistical vertex lattice 
model with weights -E^*^ (Ai, A2), is the trace of a power 
of some transfer matrix T, 



i^(Ai,A2) =Tr(T^==). 



(4) 



Here T, a (Z? — l)-dimensional tensor network itself, is 
made of all the tensors E^^ contained in some regular 
slice of the TN for F{Xi, A2), where the latter consists of 
exactly identical such slices, see Fig. ([T|). Let fia be 
the eigenvalues of T, with |/io| > > ••• > lMQ,„axl- 
Then the fidelity reads 



F(Ai,A2)= £ J2 



a=0 



SO that for large L^, and assuming |/io| > [28 



d^iXi, A2) = Mo 




(5) 



(6) 



That is, c?Ar(Ai, A2) is given by the largest eigenvalue /ig 
of T up to corrections that decay exponentially in the 
linear system size L^. Our next task is to determine /io, 
which in general depends on N, Ai, and A2. 

First we compute the left and right eigenvectors \^l) 
and \^r) of T corresponding to /io, 

(<i>i| r=($i| Mo, T = Mo |1>i?), (7) 

where we use a {D — l)-dimensional TN to represent 
them. This is achieved (again with the TN algorithms 
13, El, [3) by exploiting the fact that, e.g., ^ 
limp^oo TPl^Po) for an arbitrary state |^'o) such that 
(^o|*&_r) 7^ 0. After normalizing the states so that 
= 1, we obtain mo from 



Mo = {^L\T\<i>R), 



(8) 



by evaluating a (D - l)-dimensional TN for {^l\T\^b.), 
see Fig. ([1]). At this point, we notice that we can use 
the techniques that we have just discussed in order to 
evaluate this new TN, by reducing the calculation to a 
{D — 2)-dimensional TN, and so forth. 

We illustrate the procedure with two simple cases: (i) 
periodic chains, D = I; (ii) periodic square lattices, D = 
2. 

Case (i) [i^: Each ground state is represented as 
an MPS that consists of N copies of the tensor A^^, 



with one physical index s and two bond indices a and 
(3. The zero-dimensional transfer matrix T is given by 
i?-^(Ai,A2) = Es«'/3'(A2))*A^^(Ai), and its diago- 
nalization produces the eigenvalues MOj ■ ■ ' iMamax- 

Case (ii): Each ground state is represented as a PEPS 
on a torus with N = x Ly sites, see Fig. ([T]). The one- 
dimensional transfer matrix T is a matrix product oper- 
ator (MPO) with tensors given by the statistical weights 
E^''^ of Eq. ([3]). Its left and right eigenvectors |<I>l) 
and with maximal eigenvalue mo eh'c represented as 
MPSs with tensors L^^ and i?*^. The zero-dimensional 
transfer matrix T' reads 



tl 

£7 



y^(Z/g^)*ige'a'7's-Re"y' 



(9) 



where e = (e, e', e") and 7 are composite indices. Let Mo 
be the largest eigenvalue of T'. Then, up to corrections 
that vanish exponentially fast in and Ly, we have 



i^(Ai,A2) 



ix 
Mo 



(Mo) 



SO that 



rf7v(Ai, A2) w Mo' rf(Ai,A2) = Mo 



(10) 



(11) 



for the finite and infinite cases, respectively. 

Example: the two-dimensional quantum Ising model 
with transverse and parallel magnetic fields. As a test, 
we compute the fidelity per lattice site d(Ai,A2) for the 
two-dimensional quantum Ising model in the thermody- 
namic limit, as described by the Hamiltonian 



[r\ 



(12) 



(f,r') 



Here ax and cri ' are the Pauli matrices at the lattice site 
r, with the control parameters A and e being the trans- 
verse and parallel magnetic fields. For e = 0, the system 
has a second order phase transition at Ac « 3.044 [30| . 
whereas for A < Ac, a first order phase transition occurs 
when e changes sign. We plot d(ei,e2) and c?(Ai,A2) in 
Fig- 121 as computed from the infinite PEPS algorithm [ll| 
with bond dimension 2. We can clearly identify the first 
and second order phase transitions by a discontinuity in 
d(ei, 62) and a pinch point in d{Xi, A2), respectively. The 
two stable fixed points at A = and X — 00 are also 
characterized as the global minima of d{Xi, A2). 

Summary and outlook. The fidelity per site fi(Ai,A2) 
allows to determine the zero temperature phase diagram 
of a quantum lattice system without prior knowledge of 
order parameters. Here, we have shown how to compute 
d(Ai,A2) in the context of the TN algorithms of Refs. 

mm- 
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We envisage that this approach will become 
a preferred strategy to scan a quantum lattice system for 
possible phases and phase transitions, perhaps as a first 
step of a more comprehensive method that will subse- 
quently characterize each phase in terms of order param- 
eters, etc. An interesting question is to see whether or 



4 




FIG. 2: (color online) Fidelity per lattice site for ground states 
of the two-dimensional quantum Ising model, Eq. (|12|l , which 
is one along the diagonal. Left: for A = 2.5 (i.e. A < Ac), 
d(ei, 62) displays a discontinuity at the lines ei = and 62 = 0, 
which indicates the presence of a first order phase transition. 
Right: for e = 0, d(Ai, A2) has a pinch point at (Ac, Ac), indi- 
cating the presence of a second order phase transition. 

not the scheme works for systems with topological orders. 
On the other hand, further work is needed to perform fi- 
nite size scaling and extract the correlation length critical 
exponent by exploiting the finite TN algorithms, which 
is currently under investigation. 
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